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In  order  to  model  accurately  heat  transfer  in  PEM  fuel  cell,  a  particular  attention  had  to  be  paid  to  the 
assessment  of  heat  sources  in  the  cell.  Although  the  total  amount  of  heat  released  is  easily  computed 
from  its  voltage,  local  heat  sources  quantification  and  localization  are  not  simple.  This  paper  is  thus 
a  discussion  about  heat  sources/sinks  distribution  in  a  single  cell,  for  which  many  bold  assumptions 
are  encountered  in  the  literature.  The  heat  sources  or  sinks  under  consideration  are:  (1)  half-reactions 
entropy,  (2)  electrochemical  activation,  (3)  water  sorption/desorption  at  the  GDL/membrane  interfaces, 
(4)  Joule  effect  in  the  membrane  and  (5)  water  phase  change  in  the  GDL. 

A  detailed  thermodynamic  study  leads  to  the  conclusion  that  the  anodic  half-reaction  is  exothermic 
(A  S%v  =  -226 J  mol-1  K-1),  instead  of  being  athermic  as  supposed  in  most  of  the  thermal  studies.  As  a 
consequence,  the  cathodic  half-reaction  is  endothermic  (AS£e„  =  +62.8JmoP1  K_1),  which  results  in  a 
heat  sink  at  the  cathode  side,  proportional  to  the  current.  In  the  same  way,  depending  on  the  water 
flux  through  the  membrane,  sorption  can  create  a  large  heat  sink  at  one  electrode  and  an  equivalent 
heat  source  at  the  other.  Water  phase  change  in  the  GDL  -  condensation/evaporation  -  results  in  heat 
sources/sinks  that  should  also  be  taken  into  account.  All  these  issues  are  addressed  in  order  to  properly 
set  the  basis  of  heat  transfer  modeling  in  the  cell. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Thermal  modeling  is  of  great  interest  for  improving  hydrogen 
PEM  fuel  cell  electrical  performances.  As  highlighted  in  several 
publications  [1-3],  temperature  effects  on  fuel  cell  voltage  are  not 
negligible.  Actually,  mass  and  charge  transfers  depend  strongly  on 
temperature.  On  the  other  hand,  thermal  transfer  also  depends  on 
mass  and  charge  transfer.  Thus,  pertinent  thermal  models  consti¬ 
tute  a  key  factor  for  understanding  fuel  cell  electrical  behaviour. 

Comprehensive  reviews  of  thermal  management  issues  are  pre¬ 
sented  by  Faghri  and  Guo  [4]  and  by  Wang  [5].  Since  the  first 
description  of  coupled  heat  and  mass  transfer  proposed  by  Fuller 
and  Newman  [6],  many  similar  works  [7-9]  have  been  published. 
All  these  studies  allow  to  calculate  temperature  profiles  through 
the  membrane  electrodes  assembly  (MEA)  thickness. 

Heat  sources  distribution  and  quantification  determine  the  tem¬ 
perature  profiles  and  they  have  to  be  properly  estimated.  The 
mechanisms  responsible  for  heat  production  in  PEM  fuel  cell  MEA 
are  complex,  which  makes  it  difficult  to  localize  and  quantify  accu¬ 
rately  all  sources.  For  instance,  Ju  et  al.  [10]  and  Naterer  et  al. 


*  Corresponding  author.  Tel.:  +33  04  79  75  88  20;  fax:  +33  04  79  75  81  44. 
E-mail  address:  julien.ramousse@univ-savoie.fr  (J.  Ramousse). 


[11]  paid  a  particular  attention  to  the  estimation  of  heat  sources. 
Their  models  take  into  account  the  Joule  effect,  caused  by  the 
charge  transfer  resistances,  the  reversible  heat  of  reaction  and  the 
irreversible  heat  production  caused  by  the  electrochemical  reac¬ 
tions.  However,  the  anodic  half-reaction  is  assumed  isothermal  and 
water  sorption/desorption  phenomena  are  not  taken  into  account, 
as  in  most  of  the  models  proposed  in  the  literature.  Such  assump¬ 
tions  lead  to  inaccurate  descriptions  of  heat  transfer  in  the  cell. 
A  complete  set  of  thermodynamical  equations  was  established  by 
Kjelstrup  and  Rosjorde  [12]  who  evaluate  the  entropy  produc¬ 
tion  rate  within  the  cell.  However,  the  local  values  attributed  to 
each  phenomenon  do  not  appear  explicitly  because  of  the  global 
approach  they  used. 

This  paper  focuses  on  heat  sources  modeling  in  PEMFC.  The  next 
section  shows  how  the  total  amount  of  heat  produced  by  the  cell 
can  be  computed  from  its  output  voltage.  Then,  a  detailed  dis¬ 
cussion  about  heat  sources/sinks  quantification  and  localization 
is  given  in  order  to  properly  distribute  heat  sources  in  the  cell. 
The  following  phenomena  are  considered:  half-reactions  entropy, 
electrochemical  activation,  water  sorption/desorption  at  the  gas 
diffusion  layer  (GDL)/membrane  interfaces,  Joule  effect  in  the  mem¬ 
brane  and  water  phase  change  in  the  GDL.  Finally,  the  magnitudes 
of  the  heat  sources  are  plotted  as  functions  of  the  current  den¬ 
sity. 


0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
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Nomenclature 

a 

activity 

E 

voltage  (V) 

e 

elementary  charge  (e  =  1.6  x  10-19  C) 

E 

Faraday’s  constant  (7=  96,485  C) 

RH 

relative  humidity 

i 

current  density  (Am-2) 

Em 

membrane  thickness  (m) 

Nh2o 

water  flux  (mol  m-2  s-1 ) 

Na 

Avogadro’s  number  (NA  =  6.02  x  1023  mol-1 ) 

P 

pressure  (Pa) 

Psat 

saturation  pressure  (Pa) 

Q 

energy  (J  mol-1 ) 

Q 

surfacic  heat  source  (W  m-2) 

Q 

volumic  heat  source  (Wm-3) 

Pohm 

membrane  resistance  (Q  m2) 

reff 

effective  ionic  radius  (A) 

7 

temperature  (I<) 

W 

electrical  power  density  (W  m-2) 

Greek  letters 

AH 

enthalpy  (kj  mol-1 ) 

AG 

free  Gibbs  energy  (kj  mol-1 ) 

AS 

entropy  (J  mol-1  K-1 ) 

so 

absolute  permittivity  (F  m-1 ) 

Sr 

relative  static  permittivity  (F  m-1 ) 

X 

water  content 

°H+ 

protonic  conductivity  (S  m-1 ) 

r\ 

activation  overpotential  (V) 

Subscripts 

e~ 

electron 

H+ 

proton 

h2 

hydrogen 

h2o 

water 

02 

oxygen 

AV 

voltage  losses 

J 

Joule 

act 

activation 

reac 

reaction 

rev 

reversible 

irr 

irreversible 

sorp 

sorption/desorption 

c/e 

condensation/evaporation 

water 

cooling  water 

Superscripts 

0 

standard  conditions 

a 

anodic 

c 

cathodic 

th 

theoretical 

tot 

total 

2.  Total  heat  flux 

Heat  sources  are  directly  related  to  voltage  drops:  the  chem¬ 
ical  energy  supplied  by  the  gases  not  converted  into  electricity  is 
released  as  heat.  Fig.  1  illustrates  the  links  between  heat  production 
and  voltage  losses. 

The  chemical  power  supplied  by  the  gases  is  a  function  of  the 
enthalpy  of  the  global  reaction: 


(1) 


Polarization  curve 


Fig.  1.  Polarization  curve  and  heat  sources  (standard  conditions). 


In  standard  conditions: 

AHtot(r°,  P°)  =  Ah(,20  -  A/4  -  1/2A Hf  (2) 

From  the  reaction  enthalpy  A Hfof(7,P)),  it  is  possible  to  define  a 
reference  potential  Etot,  given  by 

Ftot  A  Htot(T,P) 

E  = - 2 F -  (3) 

In  standard  conditions,  Etot  =  1.25  V  for  vapour  water  and  Etot  =  1.48  V 
for  liquid  water  (or  1.47  V  at  P=  1  bar  and  T=  353  K).  It  represents 
the  (Active)  potential  that  the  fuel  cell  would  reach  if  the  total¬ 
ity  of  chemical  energy  were  transformed  into  electricity.  Of  course, 
according  to  the  second  law  of  thermodynamic,  only  the  Gibbs 
enthalpy  can  be  directly  transformed  into  electricity;  the  ideal  open 
circuit  voltage  (OCV)  Eth  is 

£»,_A TOP)  (4) 

Assuming  that  hydrogen,  oxygen  and  vapour  are  ideal  gases,  the 
reaction  Gibbs  enthalpy  A Gfof(7,P)  can  be  calculated  with  the  fol¬ 
lowing  expression: 


AGtot{T,P)=AG°-  f  AStot(T,  P° )  ■  dT+RT In  (  aH2°  )  (5) 

J 298  K  \aH2%  / 

where  az  =  1  in  liquid  phase;  az  =  Pl/P°  for  gases,  with  Pl  the  partial 
pressure  and  P°  the  standard  pressure.  AG°  stands  for  the  Gibbs 
enthalpy  of  the  reaction  in  standard  conditions  and  A Sf0t(7,P°) 
stands  for  the  reaction  entropy,  the  gases  being  at  standard  pres¬ 
sure.  Note  that  the  values  of  AHtot{T,P ),  A Gt0t(7,P)  and  AStot{T,P° ) 
differ  according  to  the  state  of  water  (liquid  or  vapour).  If  the  tem¬ 
perature  remains  close  to  298  K  (which  is  the  case  with  PEMFC), 
it  is  possible  to  consider  that  ASt0t(7,P°)  =  AS0  and  to  derive  the 
following  expressions  for  fuel  cell  thermodynamic  potential: 

RT  (a  H2ao/2\ 

Eth(T,  P)  =  1 .184  -  0.00023(7  -  298)  +  —  In  2  (6) 

2F  l  aH2o  / 


7^7,  P)  =  1.229-0.00085(7-298)+  ^hi(aH9a!/2)  (7) 

27  2  u2 

In  Eq.  (6),  it  is  assumed  that  water  is  produced  in  vapour  phase, 
while  it  is  produced  in  liquid  phase  in  Eq.  (7).  Of  course,  Eqs.  (6) 
and  (7)  lead  to  the  same  result  when  PH2o  =  Psat • 

The  amount  of  chemical  energy  that  can  be  converted  into  elec¬ 
tricity  is  partially  dependent  on  the  state  of  water  produced  by  the 


H2  +  l/202^  H20 
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fuel  cell  (liquid  or  gas):  the  difference  in  enthalpy  between  satu¬ 
rated  vapour  and  liquid  water  is  heat  only  but  the  thermodynamic 
open  circuit  voltage  Eth(T,P)  (6)  increases  (slightly  because  of  the 
logarithm)  if  the  vapour  partial  pressure  is  lower  than  the  satura¬ 
tion  pressure.  In  this  paper,  water  is  assumed  to  be  produced  in 
liquid  phase  in  the  hydrophilic  Nation  of  the  catalyst  layer,  which 
adsorbs  and  retains  liquid  water  [13]. 

The  real  open  circuit  voltage  is  lower  than  Eth(T,P).  According  to 
the  literature,  the  reasons  for  this  are  linked  mainly  to  [14,15]: 

•  Internal  currents  and  gas  cross-over  (hydrogen  mostly)  through 
the  membrane. 

•  Parasitic  half-reactions,  like  the  oxidation  of  platinum  and/or  oth¬ 
ers  linked  to  the  presence  of  impurities,  like  carbon  monoxide. 

The  various  voltage  drops  in  a  running  fuel  cell  and  their  effects 
on  the  polarization  curve  are  depicted  in  Fig.  1  [16].  At  low  current 
density,  voltage  losses  are  mostly  due  to  the  activation  of  electro¬ 
chemical  reactions.  Then,  the  linear  part  of  the  curve  is  commonly 
attributed  to  ohmic  losses  in  the  membrane  because  of  its  protonic 
resistance.  Finally,  at  high  current  density,  a  sharp  fall  of  the  volt¬ 
age  is  caused  by  mass  transport  limitations  (lack  of  reactive  gases 
at  the  electrodes  and/or  flooding).  The  electrical  power  lost  trough 
voltage  drops  is  irreversibly  degraded  into  heat,  Qav- 

An  additional  heat  source  associated  with  water  condensa¬ 
tion/evaporation  in  the  GDLs  has  to  be  taken  into  account,  if 
necessary.  In  practice,  low  gas  stoichiometry  leads  to  water  con¬ 
densation  in  the  gas  channels  and/or  in  the  GDLs.  The  heat  released 
by  the  water  condensation  (or  consumed  in  case  of  evaporation)  is 
not  directly  linked  to  the  cell  voltage.  Flowever,  when  too  much  liq¬ 
uid  water  accumulates  in  the  GDL,  flooding  of  the  electrode  occurs, 
resulting  in  critical  voltage  drop. 


i.  Flydrogen  oxidation  and  oxygen  reduction  are  at  the  origin 
of  a  heat  source  and  a  heat  sink  Qfeac  and  Qfeac,  respectively, 
proportional  to  the  entropy  of  each  half-reaction  and  to  the 
current  density. 

ii.  Voltage  drops  at  the  electrodes  due  to  the  activation  overpo¬ 
tentials  result  in  heat  sources  (Q£ct  and  Q£ct). 

o  And  heat  sources  Qsaorp  and  Qscorp  are  associated  with  water  sorp¬ 
tion  phenomena  at  the  membrane-electrode  interfaces. 
Actually,  these  heat  sources  are  distributed  through  the  vol¬ 
ume  of  the  electrodes.  Due  to  their  low  thickness  (about  10  p,m 
compared  to  the  membrane  or  GDL  which  have  a  thickness  about 
200  |jim),  they  are  frequently  considered  as  surface  sources.  Flow¬ 
ever,  it  is  to  be  noted  that  the  use  of  thin  membranes  like  Nation 
112  (50  |jim)  or  Gore  membranes  (18  pan)  makes  this  hypothesis 
less  tenable. 

•  Local  sources  qc/e  in  the  GDL  related  to  water  phase  change  (con¬ 
densation  or  evaporation),  if  necessary.  These  sources  can  be 
distributed  all  over  the  GDL,  depending  on  mass  and  heat  trans¬ 
fers. 

The  electrical  contact  resistances  appearing  at  the  interface 
between  the  various  MEA  layers  and  at  the  interface  between  the 
MEA  and  the  bipolar  plates  should  also  be  taken  into  account  if  they 
are  significant.  Experimental  measurements  of  the  thermal  contact 
resistances  between  Toray  carbon  paper  and  smooth  aluminium 
bronze  versus  the  compression  pressure  are  given  in  [17]. 

3.1.  Joule  effect 

Proton  transport  from  the  anode  to  the  cathode  through  the 
membrane  results  in  enthalpy  variation  due  to  the  membrane  con¬ 
ductivity: 


3.  Localization  and  quantification  of  the  heat  sources 


Qj  =  2A  Hch+  -  2  A  H“  + 


(8) 


In  this  section,  the  localization  and  the  quantification  of  each 
type  of  heat  source  (or  sink)  are  discussed.  From  a  spatial  point  of 
view,  three  cases  can  be  distinguished  (Fig.  2): 

•  Volumetric  sources  qj ,  caused  by  Joule  effect  due  to  the  protonic 
resistance  of  the  electrolyte  and  distributed  trough  the  membrane 
thickness.  The  electrical  resistance  of  the  GDLs  and  of  the  bipolar 
plates  is  often  neglected  because  of  the  high  conductivity  of  these 
carbonaceous  materials. 

•  Sources  that  are  assumed  located  at  the  electrodes: 

o  The  electrochemical  conversion  contributes  to  heat  production 
for  two  reasons: 


Anodic  Cathodic 

Bipolar  Anodic  actlvela>er  active  layer  Cathodic  Bipolar 

plate  backing  layer  Membrane  backing  layer  plate 


Oc  =QC  +QC  +  Qc 

ZZreac  act  *^sorp 


Fig.  2.  Heat  sources  distribution  in  the  cell. 


The  corresponding  voltage  loss  is  then  written: 

^ohm  —  Rohmi  —  2p  (9) 

The  proton  transfer  resistance  Rohm  is  often  estimated  after  a  pre¬ 
liminary  study  of  mass  transport  giving  water  distribution  in  the 
membrane.  Many  correlations  are  proposed  in  the  literature  for 
estimating  the  local  membrane  conductivity  from  its  temperature 
and  water  content.  The  most  used  [18-22]  are  listed  below: 


[18]: 

orH+(r)  =  (0.0013A3  +  0.0298A2  +  0.2658A.)exp 


with  Ea  =  2640  exp(-0.6A.)  + 1183 
[19]: 


orH+(T)  =  (0.5139A.  -  0.326)exp 

[20]: 

crH+(T)  =  (0.58A  -  0.5)exp 

[21]: 


1268  (303  T 


(--- 
\303  T 

(10.1) 


(10.2) 


<th+(T)  =  (0.46A  -  0.25)exp 

[22]: 


1190 


1268  (303  t)_ 

\  298.15  77  J 


<th+(T)  =  0  A.  <1.253 

erH+(T)  =  0.5738A  -  0.7192  A.  >  1.253 


(10.3) 


(10.4) 


(10.5) 
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The  total  membrane  resistance  is  then  obtained  by  integration 
of  the  above  local  conductivity  in  the  membrane  thickness: 


Rohm  ~ 


pLm 

JO 


1 

+(T,z) 


dz 


(11) 


Finally,  the  heat  source  generated  by  Joule  effect  is  computed 
from  the  membrane  resistance  f?ohm,  and  the  hydrogen  molar  flux 
Nh2  =  i/2F: 


a  =  2i  _L  =  Rohmi2 
qj  Lm  2 F  Lm 


(12) 


According  to  this  expression,  the  heat  source  generated  by  Joule 
effect  is  uniformly  distributed  in  the  membrane  thickness. 


3.2.  Heat  released  by  the  half-reactions 

Even  in  the  case  of  an  ideal  fuel  cell,  there  is  a  minimum  part  of 
chemical  energy  that  cannot  be  directly  transformed  into  work.  It 
corresponds  to: 

0th  =  AGtot(T,  P )  -  A Htot(T,  P)  =  -TAS™(T,  P)  (13) 

The  reversible  heat  released  by  the  global  reaction  of  hydrogen 
oxidation  (Eq.  (1 ))  is  a  function  of  the  reaction  entropy: 

A Sfj*  =  — 163.2 J  mol"1  K"1  at  T  =  353  K 


and  Ph2  =  Po2  =  1  bar  (14) 

Taking  gas  dilution  into  account  leads  to  the  following  expression: 

A S™(T,  P )  =  A S™(T,  P°)  -  R  ln(aH2  •  aG2)  (15) 

The  total  amount  of  heat  released  by  the  global  reaction  (in  J  mol-1 ) 
is  thus: 

Q^  =  -TASfg(T,P)  (16) 

And  the  corresponding  heat  flux  is: 

Qtot  f-m 

y<reac  ^  n  'Zreac  v 1 '  ) 


The  expression  of  Eq.  (15)  implies  that  water  is  produced  in  liquid 
phase.  If  the  water  partial  pressure  in  the  GDL  is  lower  than  the 
saturation  pressure,  water  remains  in  gas  phase  and  it  is  necessary 
to  take  account  of  the  sorption  enthalpy  at  the  cathode. 

Although  the  global  reaction  entropy  is  well  known,  determining 
which  part  has  to  be  attributed  to  each  electrode  is  a  less  simple 
task.  Indeed,  the  equations  of  thermodynamic  equilibrium  at  each 
electrode  (for  the  half-reactions  of  hydrogen  oxidation  and  oxygen 
reduction)  use  the  value  of  the  entropy  of  the  charged  species :  ASH+ 
and  ASe- .  These  values  seem  to  be  roughly  known.  Actually,  it  is 
impossible  to  make  a  solution  of  cations  without  anions,  so  that 
the  entropy  of  ions  is  not  easy  to  estimate. 

A  review  of  the  literature  highlights  this  difficulty:  as  shown  in 
Table  1,  the  values  are  widely  dispersed.  According  to  the  authors, 
the  half-reaction  of  hydrogen  oxidation  can  be  endothermic 
[25-27]  ( Q?eac  <  0)  or  exothermic  [23,24]  (Qfeac  >  0).  A  particu¬ 
lar  attention  should  be  given  to  the  conditions  in  which  these 
results  were  obtained:  on  the  one  hand,  some  authors  characterized 


Table  1 

Entropy  of  the  half-reaction  of  hydrogen  oxidation  (review). 


Entropy  of  anodic  reaction  (J  mol-1 1<-1 ) 

References 

-133.2  ±2.8 

[23] 

-42.5 

[24] 

+0.104 

[25] 

+69.1 

[26] 

+84.7  ±3.5 

[27] 

the  whole  half-reaction  of  hydrogen  oxidation  [23,27]  or  oxygen 
reduction  [25].  On  the  other  hand,  others  measured  the  entropy  of 
hydrogen  adsorption  on  Pt/C  electrodes  [26]. 

However,  Bockris  and  Conway  [28]  proposed  a  detailed  study  of 
proton  solvation  in  water  that  helps  to  estimate  the  entropy  change 
of  the  half-reaction  of  hydrogen  oxidation.  The  anodic  process  can 
be  written: 


H2gaz+2H20/Iq  2H30nq+  +  2e  (18) 

The  overall  entropy  change  can  thus  be  expressed  as 

ASarev  =  2AS“'0+  +  2ASe-  -  ASH2gaz  -  2ASh2o„,  (19) 

where  ASe-  is  the  entropy  of  the  electron  at  the  Fermi  level.  From 
Sommerfeld’s  treatment  of  the  Fermi-Dirac  statistics,  this  is  known 
to  be  equal  to  the  electronic  contribution  to  the  specific  heat  of 
metals.  It  is  about  0.25  J  mol-1  K-1 ,  which  is  far  less  than  the  degree 
of  error  in  the  other  quantities.  It  may  therefore  be  neglected. 

The  common  values  of  gaseous  dihydrogen  and  liquid 
water  entropy  are,  respectively,  ASn2gaz  =  +130. 6J  moP1 1<-1  and 

ASu2oliq  =  +69.9JmoP1  K-1  [29].  The  entropy  of  solvation  of  H30+ 
in  water  is  estimated  from  its  Gibbs  enthalpy  of  solvation,  using  the 
Born  charging  equation  [29]: 


AGh3+ 


{zefNfi  /  __1_\ 
87TS0reff  V  Sr) 


(20) 


AGh+  is  the  gain  in  free  energy  resulting  from  importing  an  elec¬ 
tronic  charge  by  integration  from  infinity  (e  =  1)  into  a  medium  of 
dielectric  constant  e.  In  our  case  z  =  1 ;  e  is  the  electronic  charge,  and 
Na  is  the  Avogadro’s  number.  Using  the  standard  thermodynamic 
relationship  (Eq.  (21))  and  the  Abegg’s  empirical  equation  for  the 
dielectric  constant  (Eqs.  (22)  and  (20))  leads  to  the  following  analyt¬ 
ical  expressions  for  the  enthalpy  and  the  entropy  of  H30+  solvation 
(Eqs.  (23)  and  (24),  respectively), 


AHh'o+  =  AGh3+-T1 


|/dAGH3Q-j 


sr  =  sr(T  =  0)  ■  exp(-TI) 


(21) 

(22) 


where  the  Abegg’s  constant  is  set  to  1  =  4.8  x  10-3  K_1  [29]. 


AHh3+ 


(ze)2NA 

87te0erreff 


(23) 


Assoi  =  (ze?NA 
h30+  87T£o£rreff 


(24) 


According  to  Conway  studies  [28],  the  enthalpy  of  H30+  solvation 
is  AH-/o+  =  -381  kj mol-1.  This  value  is  fitted  via  Eq.  (23)  to  give 

an  effective  ionic  radius  rejj  of  1.82  A.  Using  this  value  in  Eq.  (24), 
we  obtain  ASj^0+  =  +22.2 JmoP1  K-1. 

Finally,  the  overall  anodic  entropy  change  is  estimated  according 
to  Eq.  (19): 

A S%v  =  -226J  moP1  K_1  (25) 


According  to  this  analysis,  the  anodic  reaction  is  highly  exother¬ 
mic.  The  cathodic  reaction  entropy  change  is  then  deduced  from 
the  undisputable  entropy  change  of  the  overall  reaction  (Eq.  (14)): 

AScm,  =  AS“(  -  AS?e„  =  +62.8J  mol"1  K"1  (26) 


The  cathodic  process  is  therefore  endothermic. 

Because  of  the  lack  of  reliable  data,  many  authors  consider  the 
anodic  half-reaction  as  isothermal  (Qfeac  =  0)  [3,30]  and  conse¬ 
quently,  the  entropy  of  water  formation  reaction  is  fully  attributed 
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to  the  oxygen  reduction,  leading  to  an  exothermic  cathodic  half¬ 
reaction.  This  assumption,  commonly  used  in  the  literature  is  in 
contradiction  with  the  above  analysis. 

The  heat  fluxes  generated  under  reversible  conditions  by  the 
half-reactions  are  given  by 

Qieac  =  Qieac^p  (27) 

3.3.  Electrochemical  activation  of  reactions 


According  to  the  activated  complex  theory,  the  electrochemical 
reactions  are  responsible  for  a  part  of  the  energy  degradation. 

QSct  =  -TASliT,  P )  and  Q£ct  =  -TAS[rr(T,  P)  (28) 

These  irreversibilities  result  in  overpotentials  at  the  electrodes, 
depending  on  current  density,  geometrical  parameters  and  reac¬ 
tants  properties  (temperature,  pressure  and  concentration). 

r)a  =  -  and  rf  =  -  (29) 

zr  zr 

Although  the  anodic  and  cathodic  overpotentials  pa  and  ric  are  often 
calculated  from  a  Tafel  equation,  they  can  be  estimated  thanks  to 
any  model  of  mass  and  charge  transfer  in  the  electrodes  [7,31-34]. 

Knowing  the  overpotentials,  the  heat  flux  generated  by  the  acti¬ 
vation  of  the  electrochemical  reactions  is  given  by 

QSct  =  -Q£ctf:  =  riai  and  Q£ct  =  -Q£ct ^  =  rfi  (30) 

3.4.  Sorption/desorption 


Since  water  is  assumed  to  be  produced  in  liquid  phase,  sorp¬ 
tion/desorption  of  the  water  must  be  taken  into  account  at  the 
electrodes.  The  sorption  phenomenon  is  governed  by  the  balance 
between  a  liquid  phase  (adsorbed)  and  a  gaseous  phase  (desorbed) 
at  an  interface.  Water  in  fuel  cells  is  subjected  to  sorption  at  the 
interfaces  between  membrane,  where  water  is  in  liquid  phase  and 
the  GDLs,  where  water  is  mainly  vapour.  The  sorption  enthalpy 
AHsorp  is 


AHS0lp  =  AH(2oad  -  A<0gas 


(31) 


In  [35],  the  authors  observed  that  the  sorption  enthalpy  of  water 
in  Nafion  117  varies  with  the  membrane  hydration  and  that  it  is 
dependent  on  pre-treatment.  However,  no  general  trends  can  be 
drawn.  Watari  et  al.  [36]  present  a  review  of  the  values  proposed 
in  the  literature.  They  also  measured  values  decreasing  from  52  to 
45  kj  mol-1  in  the  case  of  sulfonated  polyimide  membranes.  More 
recently,  Burnett  et  al.  [37]  experimented  with  Nafion  112  and  found 
values  pretty  close  to  the  latent  heat  of  water.  All  these  results  are 
listed  in  Table  2. 

Since  most  of  the  values  encountered  in  the  literature  are  close 
to  the  latent  heat  of  water  (L„  =  41.6  kj  mol-1  at  353  K  [29]),  the 
sorption  enthalpy  of  water  is  often  assumed  equal  to  Lv. 

The  heat  released  (or  consumed)  by  this  change  of  phase  is  a 
function  of  the  water  flux  subjected  to  sorption/desorption  and  of 
the  sorption  enthalpy  AHsorp.  The  water  flux  at  the  membrane/GDLs 
interfaces  can  be  estimated  according  to  any  description  of  mass 
transport  in  the  membrane  electrode  assembly  (MEA).  This  flux  is 
perpendicular  to  the  membrane  and  counted  algebraically,  which 
makes  it  possible  to  treat  sorption  and  desorption  in  the  same  way. 


Q-sorp  —  AHsorpN}\2o  an^  Qforp  ~  —  AHsorpN^Q  (32) 

withN^Q  =  N“20  +  i/2F  to  take  into  account  of  the  water  produced 
in  liquid  phase  at  the  cathode.  Note  that  in  some  particular  cases, 
the  water  adsorption  and  desorption  in  the  membrane  can  generate 
important  heat  source  and  sink  at  the  electrodes  (for  instance  when 


Table  2 

Water  sorption  enthalpy  in  membranes  (review). 


Enthalpy  sorption  (kj  mol-1 ) 

Reference 

Comments 

46 

[38] 

For  Nafion  112 

For  aw  >  0.2 

From  60  to  43 

[37] 

For  Nafion  112 

Decreasing  with  increasing 
water  uptake 

From  52.3  to  16.7 

[39] 

For  Nafion  117 

For  aw  >  0.2 

Decreasing  with  increasing 
water  uptake 

28 

[40] 

For  Nafion  117 

For  aw  >  0.2 

From  45  to  110 

[35] 

For  Nafion  117 

A  function  of  water  uptake  and 
Nafion  pretreatement 

45 

[36] 

For  Nafion  117 

43 

[41] 

For  non-sulfonated  polyimide 
membranes 

From  52  to  45 

[36] 

For  sulfonated  polyimide 
membranes 

Decreasing  with  an  increasing 
water  uptake 

one  side  of  the  membrane  is  in  contact  with  dry  hydrogen  while 
the  other  is  in  contact  with  humidified  air).  This  phenomenon  is 
neglected  in  most  of  the  thermal  studies. 

3.5.  Water  condensation! evaporation 

Water  condensation  occurs  when  the  partial  pressure  of  vapour 
reaches  the  saturation  pressure.  This  saturation  pressure  depends 
strongly  on  temperature  as  shown  by  Antoine’s  equation  with  the 
following  coefficients: 

Psat  =  P°  exp  (l3.669  -  5°9^'23)  (33) 

In  order  to  determine  the  condensation  sites,  the  vapour  pres¬ 
sure  profiles  have  to  be  compared  to  the  fields  of  vapour  saturation 
pressure  deduced  from  temperatures.  Due  to  water  production  at 
the  cathode,  vapour  partial  pressures  are  slightly  higher  nearby  the 
membrane  and  lower  at  the  bipolar  plates/GDL  interface.  However, 
according  to  the  results  presented  in  [42],  the  vapour  pressure  gra¬ 
dients  in  the  GDL  are  lower  than  the  vapour  saturation  pressure 
gradients  resulting  from  temperature  distribution,  even  with  a  high 
GDL  thermal  conductivity.  This  means  that  condensation  will  occur 
preferentially  at  the  cold  spots  of  the  cell  (i.e.  at  the  ribs  of  the  bipo¬ 
lar  plates),  where  the  saturation  pressure  is  the  lowest,  although  the 
cathode  remains  the  place  where  the  vapour  partial  pressure  is  the 
highest.  This  observation  is  confirmed  by  the  fully  coupled  model 
of  heat  and  mass  transfers  presented  by  Weber  and  Newman  [43] 
but  it  is  in  contradiction  with  conclusions  issued  from  isothermal 
mass  transfer  models,  where  condensation  could  only  occur  at  the 
cathode.  This  remark  highlights  the  importance  of  accurate  ther¬ 
mal  and  mass  coupled  modeling  to  determine  the  condensation 
sites  in  the  cell.  To  confirm  these  remarks,  a  complete  modeling  of 
coupled  heat,  charge  and  mass  transfers  has  to  be  carried  out  in  a 
bi-dimensional  (or  three-dimensional)  geometry. 

Once  the  condensation  sites  identified,  the  associated  heat 
source  is  easily  computed  from  the  latent  heat  of  water  LV  =  (T )  at 
the  temperature  T : 

Lv{T°)  =  AH{l20i.-AH{i20mp  (34) 

At  353  K,  Lv  =  41.6  kj  mol-1  [29]. 

Let  us  denote  nvap^uq  the  water  quantity  to  be  condensed. 
nvap^iiq  is  positive  when  condensation  occurs,  and  negative  in  case 
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Heat  released  in  the  cell 


Current  density  (A.rrr2) 

Fig.  3.  Heat  sources  as  functions  of  current  density. 

of  evaporation.  The  local  heat  source  associated  to  water  conden¬ 
sation/evaporation  is  given  by 

Qc/e  =  Wvap^liqLv  (35) 

The  same  expression  is  valid  for  both  evaporation  and  condensa¬ 
tion,  since  nvap^uq  is  counted  algebraically. 

4.  Heat  sources  variation  with  the  current  density 

RHa  =  RHC  =  0.90 

Twater  =  7h2  =  T/\ir  =  353  K 

Fig.  3  depicts  the  relative  importance  of  heat  sources  (by  ref¬ 
erence  to  the  membrane  surface:  Qj  =  qjLm ),  as  functions  of  the 
current  density.  It  is  assumed  that  no  water  condensation  occurs 
within  the  fuel  cell.  The  membrane  resistance  Rohm  is  set  to 
4.4  x  10-5  £2  m2,  which  can  be  considered  as  the  worst  case,  cor¬ 
responding  for  instance  to  a  125  pum-thick  membrane  humidified 
only  by  vapour  (without  liquid  water  on  the  electrodes)  (Eq.  (11)). 
Actual  values  of  membrane  resistance  currently  reported  in  the 
literature  can  be  about  three  times  lower  [44].  The  activation  over¬ 
potentials  are  estimated  thanks  to  a  Tafel  law.  The  aim  of  Fig.  3  is 
only  to  compare  the  relative  magnitude  of  heat  sources  of  various 
origins.  A  detailed  model  of  heat  and  mass  transfer  is  needed  for 
the  proper  evaluation  of  the  local  heat  sources/sinks  in  the  cell,  as 
functions  of  local  concentration  and  water  fluxes;  this  is  out  of  the 
scope  of  this  paper.  Unfortunately,  this  representation  conceals  the 
effects  of  local  heat  sinks. 

The  heat  source  associated  with  the  global  reaction  of  hydrogen 
oxidation  is  calculated  as  follows: 

Qleac  =  Qreac  +  tteac  ~  ^  AHsolp  (36) 

with  A HSOrp  =  Lv-  Q^ac  and  Qfeac  are  the  anodic  and  cathodic 
reversible  heat  sources  related  to  the  electrochemical  half¬ 
reactions,  respectively.  The  last  term  in  the  right  hand  side  of  (Eq. 
(36))  stands  for  the  sum  of  sorption/desorption  contributions  at  the 
electrodes  (Eq.  (32)).  The  result  of  (Eq.  (36))  is  consequently  nearly 
independent  from  mass  transport  in  the  membrane.  The  activation 
overpotentials  depend  slightly  on  hydrogen  and  oxygen  concen¬ 
tration  at  the  electrodes.  Only  the  Joule  effect  in  the  membrane 
depends  significantly  on  its  water  content. 

Fig.  3  shows  that  the  activation  overpotentials  are  the  highest 
heat  sources.  The  global  reaction  of  hydrogen  oxidation  Q^c  cannot 


be  neglected:  it  is  the  lowest  heat  source  in  the  present  example  but 
it  is  expected  to  be  more  significant  than  the  Joule  effect  in  a  fuel 
cell  using  a  thinner  membrane.  These  results  are  in  good  agreement 
with  those  presented  by  Ju  et  al.  [10]. 

As  detailed  previously,  although  the  whole  reaction  results  in 
heat  production,  the  cathodic  half-reaction  is  endothermic  and 
the  anodic  one  is  exothermic.  Thereby,  this  analysis  is  in  contra¬ 
diction  with  the  assumption  commonly  found  in  the  literature 
of  an  exothermic  cathodic  reaction  (and  isothermal  anodic  half¬ 
reaction).  In  the  same  way,  water  desorption  heat  sink  is  partly 
counterbalanced  by  the  water  sorption  heat  source.  However,  one 
can  imagine  that  in  the  case  of  large  water  flux  through  the  mem¬ 
brane,  water  sorption/desorption  would  result  in  a  large  heat  sink 
at  one  electrode  and  a  large  heat  source  at  the  other  one. 

5.  Conclusions 

The  global  heat  flux  generated  by  a  fuel  cell  is  directly  dependent 
on  its  electrical  performances  and  is  therefore  easy  to  estimate. 
However,  an  accurate  review  of  all  heat  sources  within  the  MEA 
put  forward  some  difficulties  in  estimating  local  heat  sources/sinks, 
which  results  in  bold  assumptions  commonly  encountered  in  the 
literature. 

First,  because  of  a  lack  of  pertinent  thermodynamical  study 
on  both  half-reactions,  most  of  the  authors  assume  that  the 
anodic  half-reaction  is  isothermal  and  that  consequently,  the 
entropy  of  water  formation  reaction  is  fully  attributed  to 
the  oxygen  reduction,  leading  to  an  exothermic  half-reaction 
(A  =  -163. 2  J  mol-1  K_1).  The  analysis  presented  in  this  paper 
rather  suggests  that  the  oxygen  reduction  reaction  is  endother¬ 
mic  (A  Srev  =  +62. 8  J  mol-1  K-1),  whereas  the  hydrogen  oxidation 
reaction  is  highly  exothermic  ( AS{?ey  =  -226 J  mol-1 1<-1 ).  As  a  con¬ 
sequence,  thermal  gradients  in  the  different  layers  of  the  MEA 
(membrane,  GDL,  etc.)  are  probably  larger  than  routinely  thought. 

The  thermal  description  also  necessitates  to  take  account  of 
water  sorption/desorption  phenomena  at  the  electrodes:  heat  sinks 
or  sources  can  appear,  depending  on  the  direction  and  intensity  of 
the  water  flux  through  the  membrane. 

Finally,  the  water  phase  change  can  result  in  local  heat 
sources/sinks  in  the  GDL  also.  Since  the  water  saturation  pressure 
is  a  function  of  temperature  (Eq.  (33)),  water  phase  change  is  gov¬ 
erned  by  both  thermal  and  water  distributions  in  the  GDL.  A  fully 
coupled  model  of  heat  and  mass  transfers  is  thus  needed  to  identify 
the  place  of  water  condensation  and  evaporation. 

As  a  conclusion,  it  seems  that  most  of  the  thermal  transfer 
studies  tend  to  oversimplify  the  models.  The  global  heat  source  is 
generally  well  identified,  but  badly  apportioned.  As  a  consequence, 
heat  sinks  are  almost  always  neglected  and  thermal  gradients,  in 
particular  in  the  membrane,  may  be  larger  than  expected. 
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